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ABSTRACT 


The  several  equations  of  motion  governing  the 
dynamic  stability  of  a  parachute -load  system,  in  which  the 
parachute  as  well  a3  the  load  possesses  aerodynamic  drag 
and  stability  characteristics,  are  established.  The  general 
equations  are  linearized,  which  process  provides  satisfactory 
results  for  relatively  small  deflections.  A  further  sim¬ 
plification  is  accomplished  under  the  assumption  of  a 
vertical  descent.  A  numerical  example  is  used  to  illustrate 
the  application  of  the  analytical* methods. 
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SYMBOLS 


dock 


slope  of  the  canopy  normal  force 
coefficient  under  static  conditions 


slope  of  the  load  normal  force  coefficient 
under  static  conditions 


Ai,A2, A3, 
Bi,B2,B3 


constants  of  integration 

normal  force  coefficient  of  the  canopy 

normal  force  coefficient  of  the  load 

tangent  force  coefficient  of  the  canopy 

tangent  force  coefficient  of  the  load 

apparent  moment  of  inertia  of  the  entire  system, 
including  the  effect  of  the  enclosed  air  mass, 
about  its  center  of  mass  (slug-ft^) 

moment  of  inertia  of  the  canopy  and  load  about 
the  center  of  mass  of  the  system  (slug-ft2) 

total  moment  of  inertia  of  system  (slug-ft2) 

I 

dimensionless  moment  of  Inertia, 

ipSc% 

distance  between  the  center  of  pressure  of  the 
canopy  and  center  of  mass  of  the  system  (ft) 

distance  between  the  center  of  volume  of  the 
canopy  and  center  of  mass  of  the  system  (ft) 

distance  between  the  center  of  mass  of  the  load 
and  center  of  mass  of  the  system  (ft) 

distance  between  the  center  of  mass  of  the 

canopy  material  and  center  of  mass  of  the  system  (ft) 

distance  between  the  center  of  pressure  of  the 
load  and  center  of  mass  of  the  system  (ft) 

dimensionless  length,  -W 

Scy2 

apparent  mass  of  the  canopy  including  the  effect 
of  the  enclosed  air  mass  in  the  x  direction  (slug) 


▼i 


apparent  mass  of  the  load  in  the  x  direction  (slug) 

apparent  mass  of  the  canopy  including  the  effect 
of  the  enclosed  air  mass  in  the  y  direction  (slug) 

apparent  mass  of  the  load  in  the  y  direction 

mass  of  the  canopy  material  (slug) 

mass  of  the  suspended  load  (slug) 


dimensionless  mass.  ; - 5 

’  2^2 

normal  force  acting  on  the  canopy  (lb) 
normal  force  acting  on  the  load  (lb) 
projected  radius  of  the  canopy  (ft) 
characteristic  area  of  canopy  (ft2) 
characteristic  area  of  load  (ft2) 
tangent  force  acting  on  the  canopy  (lb) 
tangent  force  acting  on  the  load  (lb) 
velocity  of  the  center  of  mass  of  the  system 


velocity 

(ft/sec; 


velocity  of  the  center  of  volume  of  the  canopy 
(ft/sec) 

equilibrium  velocity  of  the  system  (ft/sec) 

velocity  of  the  center  of  mass  of  the  load 
(ft/sec) 

velocity  component  in  the  direction  of  the  system 
axis  (ft/sec; 

velocity  component  perpendicular  to  the  system 
axis  (ft/sec) 


dimensionless  velocity,  ^ 

weight  of  the  parachute  material  (lb) 

weight  of  the  load  (lb) 

angle  between  the  velocity  vector  of  the  center 
of  mass  of  the  system  and  the  canopy  axis  (radians) 

angle  of  attack  of  the  center  of  volume  of  the 
parachute  (radians) 


▼ii 


angle  of  attack  of  the  center  of  mass  of  the  load 
(radians) 


ay 

Q  angle  between  the  velocity  vector  of  the  center 

'  of  mass  of  the  system  and  the  vertical  (radians) 

area  ratio  =  S /  /S0 

0  angle  between  the  canopy  axis  and  the  vertical 

(radians) 

\A2^3  roots  of  the  frequency  equation 

p  air  density  (slug/ft^) 

Vet 

T  dimensionless  time  =  — t- 

Sc2 

to  angular  velocity  of  the  system  (radians/sec) 

Subscripts 

(  )0  initial  value  of  a  quantity  (at  T  =  0) 

Superscripts 

(  )  differentiation  with  respect  to  real  time,  t 

(  ) *  differentiation  with  respect  to  non-dimensional 

time,  *t 

(  )  dimensionless  term 
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I.  INTRODUCTION 


The  present  report  is  concerned  with  an  extension 
of  the  dynamic  stability  study  of  a  parachute -point  mass 
system  considered  in  Ref  1.  Reference  1  considered  a  simpli¬ 
fied  system  composed  of  a  point  mass,  possessing  neutral 
stability  characteristics,  and  a  statically  stable  canopy. 

One  must  realize,  however,  that  for  most  practical  applica¬ 
tions,  the  load  will  possess  a  degree  of  stability  (or 
Instability)  of  its  own.  Thus,  it  appears  advantageous  to 
attempt  an  analysis  of  such  a  parachute -load  system. 

Therefore,  this  report  presents  an  analysis  in 
which  the  load  exhibits  aerodynamic  characteristics  of  its 
own.  In  addition  to  these  characteristics,  the  motion  of 
the  system  depends  upon  the  aerodynamic  coefficients,  their 
derivatives,  and  upon  the  physical,  as  well  as  apparent  mass 
and  moment  of  inertia  of  the  parachute. 

The  general  governing  equations  are  ultimately 
simplified  for  vertical  or  near  vertical  descent. 


II.  EQUATIONS  OP  MOTION 


The  motion  of  a  parachute -load  system  Involves, 
in  general,  six  degrees  of  freedom.  In  order  to  obtain  an 
analytical  solution  with  a  reasonable  amount  of  effort,  one 
must  consider  a  simplified  physical  model.  In  view  of  these 
circumstances,  the  following  assumptions  shall  be  utilized. 

1.  Phe  entire  system  constitutes  a  rigid  body. 

2.  The  mass  and  aerodynamic  forces  of  the  suspen¬ 
sion  lines  are  neglected. 

3.  The  effects  of  the  apparent  mass  of  the  para¬ 
chute  canopy  and  that  of  the  load  act  at  the 
canopy  center  of  volume  and  load  center  of 
mass,  respectively. 

4.  The  motion  is  restricted  to  the  x-y  plane. 

This  physical  model,  with  the  acting  forces,  is 
shown  in  Pig  1. 
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>£*. 

./la 


©  CENTER  OF  PRESSURE  \ 
OF  PAR/CHUTE  CANOPY 

©  CENTER  OF  MASS  OF 
CANOPY  MATERIAL 

©  CENTER  OF  VOLUME  OF 
PAR/CHUTE  CANOPY 

®  CENTER  OF  MASS  OF  SYSTB 

<S>  CENTER  OF  PRESSURE  OF  U 

®  CENTER  OF  MASS  OF  LOAD 


Fig  1.  The  Parachute  -  Load  System 


One  may  write  the  velocity  of  the  center  of  mass 
of  the  system  in  the  canopy  fixed  reference  frame  as: 


V=VXUVyJ 


(1) 


where: 


Vx  =  velocity  of  the  center  of  mass  in  the  direction 
of  the  canopy  axis 

Vy  =  velocity  of  the  center  of  mass  perpendicular 
y  to  the  canopy  axis. 

To  determine  the  equations  of  motion  of  the  parachute -load 
system,  one  may  use  Newton's  second  law,  which  may  be  expressed 
symbolically  as: 


£F=m 


dt 


(2) 


where: 


=  sum  of  all  external  forces  acting  on  the 
system 

m  =  physical  mass  of  the  system 
d^V 


=  acceleration  of  the  center  of  mass  of  the 


system  in  an  inertial  reference  frame. 

The  absolut  /  (total)  acceleration  may  be  expressed  as  (Ref  l): 

d(,)v  d®v  ,  8 

dT=d ir*"xV  <3) 

=  acceleration  with  respect  to  reference 
frame  2 

cu  =  angular  velocity  of  reference  frame  2  with 
respect  to  reference  frame  1. 

l/i  the  present  analysis,  the  canopy  fixed  reference  frame 
/s  chosen  to  be  reference  frame  2  and  cu  may  be  expressed  as 
cu  =  ©fc.  Using  Eqns  3  and  1,  one  finds: 


where : 


d(1>V 


=  (v*  -0Vy)l  ♦  (V©VX)|. 


(4) 


Using  this  relation  in  Eqn  2,  one  obtains: 

(mp  +m/ )  [(Vx-0Vy  )i  ♦  (Vy  +0Vx)j]  =  £f  , 


(5) 


% 


where:  mc  =  mass  of  the  canopy  material 

my  =  mass  of  the  suspended  load 

The  various  external  forces  acting  on  the  system 
will  now  be  considered.  The  aerodynamic  forces,  as  shown 
in  Pig  1,  may  be  expressed  as: 


Fa=-(T,+Tc)i-(N/+Nc)j  , 


(6) 


where : 


T/ 

Tc 

N/ 


tangent  force  acting  on  the  load 
tangent  force  acting  on  the  canopy 
normal  force  acting  on  the  load 
normal  force  acting  on  the  canopy. 


One  may  express  the  gravity  forces  as: 


Fg  =  (Wy  +\A£) cos 0 1  -  (W/+Wc)sin©I 


(7) 


The  apparent  mass  forces  can  be  expressed  as: 


“am 


=  -  nrtaxc^xji  “  [may^  3y/  +  mayc  SyJ  j  , 


(8) 


where :  m^ 

may/ 

maxc 


%yc 


apparent  mass  of  the  load  in  the  x  direction 

apparent  mass  of  the  load  in  the  y  direction 

apparent  mass  of  the  canopy  including  the 
effect  of  the  enclosed  air  mass  in  the 
x  direction  (Ref  2) 

apparent  mass  of  the  canopy  including  the 
effect  of  the  enclosed  air  mass  in  the 
y  direction  (Ref  2) 

acceleration  of  the  center  of  mass  of  the 
load  in  the  x  direction 


av  =  acceleration  of  the  center  of  mass  of  the 
load  in  the  y  direction 

ax  =  acceleration  of  the  center  of  volume  of 
c  the  canopy  in  the  x  direction 

ay  =  acceleration  of  the  center  of  volume  of 
c  the  canopy  in  the  y  direction. 
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The  velocity  of  the  center  of  mass  of  the  load  and  center  of 
volume  of  the  canopy  consists  of  the  velocity  of  the  center 
of  mass  of  the  system  $  plus  a  rotational  velocity  about  the 
center  of  mass  of  the  parachute -load  system.  Thus,  one  may 
write : 


/ 


%  =  V+£x  L  3f  =  Vxf  +  (Vy  ♦  L30)J 
Vc  =  9+5x(-L2i)  =  Vxi  +  (ty-L20)  J 


The  acceleration  of  the  center  of  mass  of  the  load  and  center 
of  volume  of  the  canopy  may  be  written  as : 

-  <% 
a/  =  3T 
- 

*  *  dt 

and,  using  relations  3  and  9,  the 


a,  =  (v*  -evy  -  L3©2]  1  ♦  [ty+0Vx 
lc  =  [Vx  -  0Vy  +  L20^|  1  +  (Vy  ♦  0VX 


above  equations  become: 


♦L3©]] 

-L2e]l 


(10) 


Or,  in  scalar  form,  one  may  write  the  corresponding  acceleration 
components  as : 

a*,  =V*  -evy  -L392 
3^  =\^  +  0VX  +  L30 
aXc=^x  ~©Vy  ♦^■2©2 
ay«  =  Vy  +©Vk  -L2© 

and  the  apparent  mass  forces,  in  accordance  with  Ean  8,  become: 


Fam  =  -  [maX/(V*  ~©Vy  -L302)  ♦  maxc(Vx  ~©Vy 

-  [may/(Vy  +0VX>L30)  ♦  ma*(Vy  +0VX  -L20)j  ] 
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(11) 


! 


* 


i 

! 


i 


» 
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Utilizing  relations  6,  7  and  11  in  Eqn  5  yields  two  scalar 
equations  of  motion  in  the  x  and  y  directions,  respectively: 


(mc  +rry  +rriax^+niaxc )  (Vx  -  0Vy )  = 

-  (T>  +TC )  +  (W/  +WC)  cos0  ♦  (maX/  L3  -maXcL2)e2 


(12) 


•  • 

(me  +  rry  ♦  may^ +mayc  )(Vy  +6\^)  = 

- N|  -Nc-(W!+Wc)sin0  +  (ma*L2-may4L3)0 


(13) 


A  third  equation  may  be  written  which  governs  the  rotational 
motion  of  the  entire  system.  This  equation,  the  angular  mo¬ 
mentum  equation,  states  that  the  sum  of  the  external  moments 
acting  about  the  center  of  mass  of  the  system  equals  the  time 
rate  of  change  of  the  angular  momentum  of  the  system  about 
its  center  of  mass.  This  statement  may  be  expressed  as: 

(Icj*US«  NcU-HLg  (14) 

where:  Ic  g  =  moment  of  inertia  of  the  canopy  and  load 

about  the  center  of  mass  of  the  system 

Ia  a  apparent  moment  of  inertia  of  the  entire- 

system  including  the  effect  of  the  enclosed 
air  mars  about  its  center  of  mass 

IrpLtj  =  lengths  (see  Fig  l). 

The  various  aerodynamic  forces  are  conventionally  expressed  as 

Vq^pvfs,  t,  iC^pvfs,  '15) 

In  addition  to  these  definitions,  it  is  convenient  to  make 
the  equations  of  motion  dimensionless  in  the  following 
manner: 


6 


m  = 


m 


P  sp 


L  = 


L 


Sj l12 


v  = 


va 


1= 


TP  St 


5/2 


(16) 


t  = 


Sj/2 


Ve  =  equilibrium  velocity. 

Introducing  these  relations,  along  with  Eqn  15,  into  Eqns  12 
through  14,  one  obtains  after  some  algebraic  manipulations: 


( 


mc  *ml  ■*rrV<’rnaxt)^x  =  -  CrVo2 

+  WtfWzl  cos  0 
ipV|_Sc 

+  ^  maX|^"3“  maxJ~  2  ^ 


(12a) 


9X<>  *  -cNt^|  -c$2 

7  pvl  Sc 

+(r%X2  -frVE3)e" 


(13a) 


dcg^je"  =-cN/?|-L5*  CncVc2Li  , 


(14a) 


where  the  prime  (')  indicates  differentiation  with  respect 
to  the  dimensionless  time  t  . 

The  equilibrium  velocity  is  defined  by: 

vv If  *WC  -(CTfS„  ♦CT.ScJjpVl  , 


which  may  be  written  as: 


7 


*  ?  «.'f*i**  Tt^As**#*#*^#****^**^********?^#*^^  ' 


\Nj,  +WC 
2 -P  ^  Sc 


(17) 


To  abbreviate  the  form  of  Eqns  12a  through  l4a,  let  us 
write: 

mx  =  mc  *mj  +ffiaX/+rriaxc 

my  =mc  *ff\f  +may/  +ma>i 

^TOT.  =^cg  +^a 

r  -¥ 

oc 

Ax  =f^ax^l-3  -  FfiaxcL.2 
Ay  s  rnayo  L2.  -  ma^  l.  3 


(18) 


Substituting  the  definitions  17  and  18  into  Eqns  12a  through 
l4a  yields : 


mx(v;  -^vy)=-rcT/v/  -cTc  ^(c^+rcT#)cose+Axe'2  (12b) 

my(v'  +  6\)=-tcH/tf-  cNcv?-(cTc^cT/)sine+Aye"  (13b) 

Itot©  =  ~^CN/  ty2  L5  +  C  Nc Vc2  L1  .  ( 14b ) 


In  order  to  obtain  £he  final  equations  of  motion  in  a  con¬ 
venient  form,  the  following  relationships  may  be  deduced 


8 


tan  oc t  =  tanoc 


L3Q 

Vcosoc 


tana.  =  tana  +  — 

c  Vccsa 


(23) 


Utilizing  relations  21  and  22,  the  equations  of  motion  12b 
through  14b  assume  the  following  form: 


mjv'cosa  -Va'sin cc  +  V©'sin  a  j  = 

-JCT|[v'2-2C3V©'sina  *  l^©2] 
-C  [Vi  2L;p©sincc  *  H2^] 


♦ (CT'JC  )cos©  *  A©2 

L  i*  X 


rHyJ-V'sina  -Voc'cosa  +  tf©  cos  a  ]  = 

-jrCN[7-2L3^sina+t^] 

-cJV2*2L2tf€f  sina  +  C^'2] 
-(C+fC *  )sin©  ♦  A©* 

y  1 1  y 


(24) 


(25) 


Jr©;--^C5[V?2^©'s,na  .^] 
.C^^fV2.  2L2V©'sha  *  C.2©2] 


(26) 
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III.  A  FIRST  ORDER  METHOD  OF  OBTAINING 
NUMERICAL  SOLUTIONS 


A  simplified  numerical  method  of  solution  will 
now  be  outlined.  In  the  most  general  situation,  certain 
initial  conditions  must  be  given.  A  typical  set  of  Initial 
conditons  is  given  by: 

at  x- o,  ©=e0,  a=oc0,  q'=Qo,  v=v0 


Let  us  first  consider  Eqn  23.  One  observes  that  these  initial 
conditons  are  sufficient  to  determine  the  initial  angle  of 
attack  of  the  load  and  canopy,  CC.  and  CCc  .  It  may  be 
assumed  that  the  aerodynamic  coerficients°CNc,  Ctc,  Cn*  and 
Cij,  are  known  functions  of  0Cc  and  OCj  respectively.  Con¬ 
sequently,  OC^  and  CCCq  determine  the  initial  values  of 
the  coefficients.  Substituting  the  initial  coefficient 
values  and  the  initial  conditions  from  above  into  Eqns  24 
through  2 6  yields  three  equations  of  the  form: 


Ay0'+  B,cC=Q 

(24a) 

Ay0/  +  B2a;=C2+  D& 

(25a) 

A^9o  =  C3  . 

(26a) 

The  constants  A^,  Ag,  .  .  *  are  determined  from  Eqns  24 
through  26.  For  example,  by  comparison,  one  finds  that 
At  =  %  cos  OC0 ,  A2  =  -  Wx  vo  sin  CXo t  etc.  Thus  Eqns  24a 
through  26a  represent  three  linear  equations  from  which  one 
can  determine  v0',  CC0 * ,  and  ©0". 

One  next  selects  a  small  dimensionless  time  Interval 
and  calculates: 

d?  = 

4©'  =  ©b<;r 
JOL-Q&SZ 
✓©  =  00  it 
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-rivs#  -r***  -  •W'M.w 


Thus,  after  the  time  interval  (ft  ,  the  new  values  of  the 
variables  become: 


V1  = 

vG+  /V 

e;  = 

e;w07 

<*1  = 

oc0+  <f  oc 

01  = 

CD 

One  now  has  a  new  set  of  conditions  Vi,  ©i*,OC  i,  01,  and 
the_j>re  ceding  steps  can  be  repeated  to  determine  new  values 
of  V,  ©’,  @  ,  andoc  •  Tills  procedure  may  be  repeated 
indefinitely  until  the  entire  trajectory  is  determined  as: 

V  =  1ft) 

e  =  f2(t) 
oc  =  f3(C) 

The  above  procedure  is  equivalent  to  expanding  V,  0  ,  and  oc 
in  a  Taylor  series  and  retaining  only  the  first  two  terms. 
In  this  manner,  the  nonlinear  differential  equations,  24 
through  26,  including  a  nonlinear  C^,  Ct,  -cc  relationship, 
can  be  solved. 


IV.  LINEARIZED  THEORY 


If  one  wishes  to  consider  only  the  class  of  motions 
where  the  oscillations  are  small  (  oc  small),  the  following 
equations  are  applicable: 

SIN  a  =  a  COS  cc  =  1 

In  addition,  if  the  trajectory  is  almost  vertical  ( small), 
one  may  write : 


S!Ny0  =  COS  fi  =1 

Utilizing  these  assumptions  in  Eqns  24  through  26  and 
neglecting  second  order  terms  such  as  OCOC*  ,  one  obtains: 

mj'--  -(*-CVCtHV2-1>  (27) 
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myV(e'-a')  *  -( fcNl*c^>vMcT,,rcT|>e.  Ay6- 


^TDT®  =  ^Nc*~l  "^NJl*-5^  ^ 


Equation  27  may  be  integrated  directly  to  give 
the  variation  of  V  with  -j-  . 


where  V0  is  the  value  of  V  at  x  -  0* 

Equations  28  and  29  can  be  presented  in  a  .more 
explicit  form  by  realizing  that  for  small  oscillations,  C<rt 
and  Ctc  are  constant  while  Cjj<i  and  Cnc  are  linear  functions 
of  oc^  and  oc  c  respectively.  Thue^  one  may  write: 


CNjl=  a*  a*  (: 

~  9c  CCc 

where:  a»  =  slope  of  versus  oc.  under  static 

conditions  " 

a„  -  slope  of  Cm  versus  oc  «  under  static 
conditions.0 

Introducing  these  relations  into  Eqns  28  and  29 
and  assuming  that  the  system  is  descending  at  approximately 
its  equilibrium  speed  (l.e.,  V  =  l),  one  finds: 


Frye1- a'y-&x<xr  (qc4rcT,)e  ♦  Aye^occ 


(28a) 


Itot6'  =  'acacL'  -  fa.a.L.) 


(29a) 


Again  utilizing  the  assumption  of  small  oscillations,  Eqn  23 
assumes  the  form: 


oc*  =  oc  -  tp 
occ  =  oc  +  l20/ 


(23a) 


Introducing  these  relations  into  Eqns  28a  and  29a  yields, 
after  rearranging: 


Ay©  +  [*L3a2-  L2ac-  mvJ©  -  [C-^  +  t C-J  ® 
+  myCC'-  [jra^+aJOC  =  0 


(28b) 


IT0]e"-  [yC3L5aJl+  L^ac]©'-  (L^-sO^OC  =0 


(29b) 


Equations  28b  and  29b  are  coupled  but  linear,  and 
one  may  assume  a  solution  of  the  form: 


a  =  AexX 
0  =  Be** 


(32) 


Substituting  these  functions  into  Eqn  28b  and  29b  yields, 
after  rearranging: 


where: 


a„A  +  a12B  =0 
a21  A  +  a22B  -  0 

an  =  iTiyA-fya,  +  a^j 

a12  ”  AyA  +  £  +  ^3al  *  myJ  £0rc  + 

a21  =  yLsar^iac 

a22  =  ^TOT^  "  (-3^-5ax  +  ^-2ac|  ^ 


(33) 


(33a) 
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A  nontrivial  solution  of  Eqn  33  for  A  and  B  exists 
if  and  only  if  the  determinant  of  the  coefficients  is 
identically  zero.  That  is,  if: 


all 

a12 

a21 

a22 

Expanding  the  determinant,  one  finds  after  rearranging: 


aA3  +  bA2  +  cA  +d  =0 


(34) 


The  coefficients,  a,  b,  c,  and  d,  obtained  by  a  direct 
expansion  of  the  determinant,  can  be  written  with  the  aid 
of  relation  18  and  33a  as: 

a  z^Jtot 

b  =  -  ( mc+ rn, )( + rl^L5a, ) 

-( L2+n3)(ffiay(L1  ac+  rr^Lja/)  -T-^ac+^ai) 

(35) 

C  =  #'(£|+  L-sXL^+t^SjSc  ~  rfly(Hiac  ”^-53/) 

d  =  -(Ctc-*iCt|X£i  ac-  *C5ar) 


Equation  34  is  referred  to  as  the  frequency  equation 
of  the  system.  In  general,  it  yields  three  distinct  values 
of  A  . 

Routh's  criteria  (Ref  3)  requires  that  for  a 
dynamically  stable  system,  the  following  inequalities  be 
satisfied: 

d  >0 

bod  (36) 


In  essence,  the  angle  of  attack  OC  and  the  related 
angle  ©  (Eqn  32)  decay  with  time  if  the  real  part  of  the 
roots  of  the  frequency  equation  (34)  are  negative.  This  is 
the  case  if  Routh's  criteria  is  satisfied. 


a  >  0 
b  >  0 
c  >  0 
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Now,  details  of  the  solution  of  the  governing 
equations  shall  be  discussed.  Since  three  roots  of  the 
frequency  equation  (  A]_,  A  g*  A3)  can  be  found,  the  general 
solution  of  the  problem  may  be  written  as: 


a  =  +  A^+  A3e^ 

0  =  B^  +  Bg^+Bs^ 


where  A]_,  A2,  Ao,  B^,  Bg,  and  B3  are  the  constants  of  in¬ 
tegration.  Three  relations  for  these  constants  can  be  found 
through  the  use  of  the  initial  conditions  expressed  as 
CX=0Co,  0  =©0  ,  and  0’  =  ©0’  at  X-  0.  Applications  of 
these  conditions  yield: 


A 1  +  A  2 +  A3  =  oc0 

Bi  +  B2+  B3=  0o  (38) 

A^  A2B2+A3B3=  ©o 


Additional  equations  can  be  obtained  from  either  of  Eqns 
33.  Utilizing  the  second  of  these  equations,  one  finds, 
with  the  use  of  33a,  the  three  relations: 


Aj.  .  IioTA,-(ff  L3L5aj;+  L1  Lgac  )Aj 

Bj  -  Hiac+  if  [©a^ 

where  1  assumes  the  values  of  1,  2,  and  3.  Thus,  using 
Eqns  38  and  39,  all  of  the  constants  are  specified  and  the 
solution  established, 

A  particular  numerical  example  of  a  system  com¬ 
prised  of  a  stable  parachute  and  uns table  load  is  treated 
in  Section  7. 
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(44) 


A,  -fd*-/1*#2') 

Aa-fd-V^^) 


The  general  solution  of  Eqn  40  may  now  be  written 
in  one  of  several  ways,  depending  on  the  relative  values 
and  signs  of  m  and  rJ2.  The  case  of  a  statically  unstable 
parachute,  where  ac  >  0,  is  certainly  not  capable  of  sta¬ 
bilizing  an  unstable  load.  Therefore,  this  case  may  be 
disregarded. 

Thus,  for  all  physically  realistic  systems,  one 
sees  from  cheir  definitions  that  Li,  L2,  L 3,  If  ,  and  I^qt 
are  always  positive  and  ac  is  negative.  Also,  for  most 
practical  cases,  the  slope  of  the  normal  force  coefficient 
versus  the  angle  of  attack  for  the  load  is  negative.  However, 
L5  can  be  either  positive  or  negative,  depending  on  the 
relative  sizes  of  the  parachute  and  load  as  well  as  on  the 
stability  of  the  load  by  itself. 

The  above  considerations  thus  show  that  both  m 
and  n2  can  be  either  positive  or  negative  quantities. 

However,  for  most  practical  systems,  both  m  and  n2  will 
be  positive  and  as  a  result,  only  this  case  will  be  con¬ 
sidered  in  detail  in  this  report. 

In  view  of  Eqn  44,  it  is  advantageous  to  consider 
the  following  three  cases  of  the  term  (n/m)2: 

0  <  (&f<  1 

For  this  case,  examination  of  Eqn  44  shows  that 
both  A  1  and  A  2  are  negative  and  the  general  solution 
can  be  written  as: 

0  =  A,  ex,t  +  Age*** 


Since  7\\  and  A2  are  negative,  the  angle  (©)  decays  with 
time  and  the  system  is  dynamically  stable. 


In  this  case  one  finds  that  A]_ 
giving  the  solution: 


A  2  -  -  §>  thus 
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-  ~ig»yvw  y •ysjrs,  *-f**xs^xw^!*y,-^ 


nvt 

e  =  e*  2  (A,  + Agt) 


Once  again  this  solution  indicates  a  dynamically  stable 
system. 

(m>2>1 

In  this  case,  the  two  roots  A^  and  A.  2  are  complex 
numbers  and  one  can  write: 

x  -  G1  j  00  .^~n  £  'T 

"  2  "  2  1(W 


v-**«w 

and  tbe  general  solution  can  be  written  as: 


9  =  e!?[A,sin»1/(ft)2-l‘  -t  *  Ajcosm  -^~T t 
-  “  * 


which  represents  a  stable  system.  This  case  represents  the 
most  common  situation  for  a  stable  system  and  one  can  deduce 
certain  characteristics  of  the  motion  of  the  system  from 
this  relation. 

m+ 

-  -grt 

The  damping  factor,  e  ,  indicates  the  rate  at 

which  the  oscillations  are  damped.  Large  values  of  m  correspond 
to  a  rapid  damping.  One  observes  from  Eqn  4l  that  large 
values  of  m  (rapid  damping)  correspond  to  a  combination  of 
a  very  stable  parachute  (large  negative  ac),  long  suspension 
lines,  and  a  small  moment  of  inertia,  I. 

In  order  to  solve  Eqn  46,  one  must  determine  the 
constants  of  integration,  preferably  from  the  initial  con¬ 
ditions  ©  =  ©  o,  ©'  =  ©o'  at  X  =  0.  One  finds  the  con¬ 
stants  to  be: 


1  Q> 

6'  = 

A  - 

eU?' 

Wf 

a2 

=  9. 
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1  »  f 


and  the  solution  of  Eqn  46  amounts  to: 


e  =^r 


*©<.cos^ 


-1 


-1 


(47) 


VI.  CONCLUSIONS 


With  the  assumption  of  small  oscillations  and 
near  vertical  descent,  one  may  completely  specify  the  motion 
of  a  parachute -load  system  in  which  the  load  possesses 
distinct  aerodynamic  stability  properties.  This  rigorous 
solution  requires  lengthy  calculations.  To  alleviate  this 
problem,  one  may  choose  to  solve  the  simplified  equation 
presented  in  Section  5.  This  approach  is  justified  for  f 

vertical  descent  { j3  =  0).  1 


A  substantial  portion  of  the  calculations  can  be 
eliminated  if  one  merely  wishes  to  determine  whether  or  not 
the  system  is  dynamically  stable.  The  answer  to  this  question 
is  furnished  by  Routh’s  criteria  (Section  4).  One  must. 

In  this  case,  only  determine  the  value  of  the  coefficients 
of  the  frequency  equation  (Eqns  34  and  35)  and  check  to  see 
if  they  satisfy  relations  36.  Because  of  the  number  of 
parameters  involved  and  the  way  in  which  they  enter  the 
frequency  equation.  It  is  difficult  to  draw  general  con¬ 
clusions  pertaining  to  the  stability  of  a  general  system. 
Therefore,  for  any  practical  case,  Routh’s  criteria  must  be 
satisfied  in  its  various  aspects. 

It  also  can  be  seen  that  a  system  consisting  of 
a  stable  or  unstable  load  combined  with  a  stable  parachuoe 
can  be  analyzed  in  the  same  manner  as  a  system  with  a  point 
load  and  a  stable  parachute (Ref  2).  For  the  case  of  the 
point  load,  the  governing  equations  of  this  study  may  be 
reduced  to  the  same  form  as  given  in  Ref  2  when  the  aero¬ 
dynamic  coefficients  of  the.  load  approach  zero. 
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VII.  NUMERICAL  EXAMPLE 


A  numerical  example,  concerning  a  parachute -load 
system  (Pig  2)  which  consists  of  a  five-foot  ribless  guide 
surface  parachute  canopy,  having  a  nominal  porosity,  of 
70  ft3/ft2-min,  and  a  one-foot  diameter  ogive  cylinder 
weighing  350  lbs  shall  be  used  to  illustrate  the  presented 
theory. 


The  mass  of  the  parachute,  mc,  shall  be  the  same 
as  in  Ref  2,  namely: 

mc  =  2.82  ti  r2x  10'3  slugs 

which,  in  dimensionless  terms  is: 


mc  =  0.537 


Prom  the  same  source,  the  mass  of  the  enclosed  air  amounts 
to : 


rrij  =  0.419 

and  in  dimensionless  form: 


m  j  =  0.468 

Prom  Ref  4  one  finds: 


15*?  =  0.3 
rrij 


and  therefore: 


nnaxc  =  0.1404 


The  dimensionless  mass  of  the  load  is: 

m2  =  105.17 

From  the  geometry  and  various  physical  masses  of  the  com¬ 
ponents  of  the  system,  one  finds: 


L,  =  2.2166 


L2  =  2.1665 


/ 


21 


rni‘KK'*n*v*~  —■ 


L3  =0.0064 
U4  =2.2105 


A^.c-ming  that  the  distribution  of  mass  in  the  ogive  is 
uniform,  the  moment  of  inertia  amounts  to: 

\ 

Tcg =  10.8536 

Again,  from  Ref  2  one  finds  the  apparent  moment  of  inertia 
of  the  canopy  about  the  center  of  mass  of  the  system  to  be: 


Ia  =  0.1874nr3pL| 


and  for  this  particular  configuration: 

Ta  =13204 

thus  giving  a  total  inertia  of: 


W  =12.1740 


The  center  of  pressure  of  the  ogive  cylinder  has 
been  experimentally  determined  through  wind  tunnel  tests 
and  it  is  found  to  be  located  a  distance  of  0.3 Jl  behind 
the  tip  of  the  ogive  {Note:  i  represents  the  length  of 
the  ogive  cylinder) ,  With  the  preceding  geometry,  one  can 
determine  the  distance  L 5,  in  dimensionless  form,  as: 

L5  =  0.2095 


Also,  the  above-mentioned  experiments  have  shown  that: 

CT,  =0.23 

at  =-2.55  per  radian 


and  from  Ref  2  for  the  parachute: 

CTc  =1.08 

ac  =  -0.676  per  radian 
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as : 


All  of  the  preceding  results  may  be  summarized 


mc  =0.537 
maxc=mayc=  0.1404 
m,=  105.17 
U=  2.2166 
L2=  2.1665 
C3=  0.0064 
L4=  2.2105 
L5  =  0.2095 
Ttot=12.1740 
Cy^O.23 

CTc  =  1.08 

aa  =  -  2.55  per  radian 
ac  =  -0.676  per  radian 
V  =  0.04 


In  addition  to  these  values,  it  will  here  be 
assumed  that  the  apparent  mass  of  the  load  is  negligible 
and  m8LXt  =  may£  =  0. 

Utilizing  the  above  values  in  the  frequency  equa¬ 
tion  (Eqn  34),  one  finds: 

1285  A3  +  352.47  A2  +  156.7314  A  + 1.6088  =0 


of,  dividing  by  the  coefficient  of  A^: 


A3  ♦  0.2742  A2  +  0.1219  A  +  .001251  =  0 


(48) 


Solving  this  relation  by  the  method  presented  in  Ref  5, 
one  finds: 


* 
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\  =  -0.010535 

A2  =  -  0.1236  +  0.315554  i 

A3  = -0.1286  -  0.315554 i 

(49) 


As  initial  conditions,  we  choose  at  't  -  0: 

Go  -10°=  0.1745 

©o=0 

cq,  =10°  =  0.1745 

Using  Relations  3$  and  39>  one  finds  the  values  of  the 
constants  A-x,  A^,  A3,  B^,  Bg,  B^  to  he: 

At  =-0.00008294 
A2=  0.09460  e0-39561' 

A3=  0.09460  ea39561i 

B,  =  -0.0037315  (50) 

B2=  0.09621  ea38630i 
B3  =  0.09621  e038630' 

After  several  algebraic  manipulations,  the  linearized  general 
theory  provides  the  angles  0  and  OC  as: 

©  = -0.003731  e'Q01054T 

+Ot19242  e'0-12867  cos(0.3156  7  -  0.38630)  (51) 

a  =  - 0.00008294  e  00,054  r 

+  0.1892  <s°.i286*  cos  (0.3156  7  -  0.39561) 


(52) 


Using  the  simplified  analysis  for  the  vertical  descent 
{ jB  =  0),  Eqn  47  gives: 


©  =  0.1885  e°'1304*cos( 0.3187 1  -  0.38834)  (53) 


The  relations  51  and  53  are  graphically  presented 
in  Pig  3  and  tabulated  in  Table  1.  It  can  be  seen  that  both 
approaches  provide  nearly  identical  results  and  it  appears 
that,  for  many  practical  cases,  the  simplified  method  would 
be  entirely  satisfactory.  However,  for  a  case  involving  a 
trajectory  with  a  strongly  changing  inclination  angle  or  a 
system  with  nonlinear  aerodynamic  coefficients,  a  numerical 
solution  of  the  nonlinearized  equations  may  be  required. 
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Fig  3  0  as  a  Function  of  Time  for  the  Exact  and 

Approximate  Solutions 
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13  ABSTRACT 


The  several  equations  of  motion  governing  the  dynamic  stability  of  a 
parachute- load  system,  in  which  the  parachute  as  well  as  the  load 
possesses  aerodynamic  drag  and  stability  characteristics,  are  established. 
The  general  equations  are  linearized,  which  process  provides  satisfactory 
results  for  relatively  small  deflections.  A  further  simplification  is 
accomplished  under  the  assumption  of  a  vertical  descent,  A  numerical 
example  is  used  to  illustrate  the  application  of  the  analytical  methods. 
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